Current-controlled magnetization dynamics in the spin-flip transistor 
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The current driven magnetization dynamics of a thin-film, three magnetic terminal device (spin- 
flip transistor) is investigated theoretically. We consider a magnetization configuration in which 
all magnetizations are in the device plane, with source-drain magnetizations chosen fixed and an- 
tiparallel, whereas the third contact magnetization is allowed to move in a weak anisotropy field 
that guarantees thermal stability of the equilibrium structure at room temperature. We analyze the 
magnetization dynamics of the free layer under a dc source-drain bias current within the macrospin 
model and magneto-electronic circuit theory. A new tunable two-state behavior of the magnetization 
is found and the advantages of this phenomenon and potential applications are discussed. 

I. INTRODUCTION 

The current induced magnetization excitation predicted by Slonczewski and Berger 1,2 has attracted considerable 
attention and the prediction of current-induced magnetization reversal has been confirmed by many experiments 
in nano-pillar structure consisting of two ferromagnetic layers with a high ("fixed") and a low ("free") coercivity, 
separated by a normal metal spacer^^ii Meanwhile, the investigations of charge and spin transport in thin-film 
metallic conductors structured on a planar substrate have also been carried ou tiL£iiLi2iiiii£ The advantages of the 
planar structures are the flexible design and the relative ease to fabricate multi-terminal structures^ Recently, non- 
local magnetization switching in a lateral spin valve structure has been demonstrated^* In the present article we 
present a theoretical study on the dynamics of a lateral spin valve consisting of a normal metal film that is contacted 
by two magnetically hard ferromagnets. As sketched in Fig. ^ a slightly elliptic and magnetically soft ferromagnetic 
film is assumed deposited on top of the normal metal to form a spin-flip transistor^ The magnetization direction of 
the source-drain contacts lies antiparallel to each other in the plane of the magnetization of the third (free) layer. The 
configuration in which the source-drain contact magnetizations are oriented perpendicular to the plane is considered 
elsewhere^ A convenient and accurate tool to study the dynamic properties of our device is the magnetoelectronic 
circuit theory (MECT) for charge and spin transport^ combined with the Landau-Lifshitz-Gilbert equation in the 
macrospin model. The spin flip scattering in normal and ferromagnetic metals and the spin-pumping effect are also 
taken into account JLii 

The article is organized as the follows: In fjnj we briefly review the MECT and Landau-Lifshitz-Gilbert equation 
for the macrospin model. In flIIII calculations of the spin transfer torque for our device are presented, ill VI is devoted 
to the discussion of thermal (in)stability and in JV] the magnetization dynamics is treated. The conclusions are 
summarized in WII 
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FIG. 1: The model system contains a normal metal sandwiched by two ferromagnetic leads and a circular soft ferromagnet 
film (e.g., permalloy) on top of the normal metal. The magnetization unit vectors mi, ni2, and 1H3 are initially aligned in the 
same, i.e., x — y plane. 
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II. MAGNETO-ELECTRONIC CIRCUIT THEORY 



We first consider a ferromagnet-normal metal (F|N) interface at quasi-equilibrium. The ferromagnet at a chemical 
potential [j,q and spin accumulation fifm aligned along the magnetization direction. The chemical potential and spin 
accumulation in the normal metal are denoted by /ig and vector S. The charge current I c (in the unit of Ampere) 
and the spin current I s (in the unit of Joule) entering the normal metal node are given by^ 

Ic [2s(Mq -l*o)+ P9Vs - P9 (S • m)] 
I, =f[2p04 U$) + /if (1 - »Jr)(S • m)]m 

-M S -M<s> m >. (1 ) 

From eq. we may project out the component of I s that is perpendicular to the magnetization direction and governs 
the spin transfer torqu e)^ 15 

mxl s xm=^[S-(S-m)m] + ^(Sxm) . (2) 

07T 07T 

In the above notations, the dimensionless total conductance g = g' + g^ and the mixing conductance are given by 
Landauer-Biittiker formulae, i.e., 

„TU) _ m - V \r nm 2 

nm 

g n = m -j2 r v n ( r i m y ■ ( 3 ) 

nm 

where rKJ is the probability of a spin up(down) electron in mode m reflected into mode n in the normal metal and 

M is the total number of channels. The contact polarization is defined by p — (g" — g^)/{g^ + g^)- The pumping 
current generated by the motion of the magnetization isi£ 

■ (4) 



8tt V " dt ' dt , 

We consider the situation in which the dimension of the normal metal is smaller than the spin-flip length, so that the 
spin accumulation does not vary spatially in the node. 

!(/) = £S (5) 

where <?/ = hi/£>os^ N /tj ' , ^dos an d Vat are the density of states of the electrons and the volume of the normal 
metal, is the spin-flip relaxation time in the normal metal node. The charge and spin currents entering the normal 
metal obey the conservation laws 

£/ C)i = o, 

i 

•■:::) (6) 



III. SPIN TRANSFER TORQUE 

In the structure depicted in Fig. ^ the source-drain magnetizations are aligned anti-parallel along the y-axis in order 
to inject a large spin accumulation into N, i.e., mi = (0, +1, 0) and ni2 = (0,-1, 0). Connecting the ferromagnets to 
reservoirs and applying a bias current Iq via the two ferromagnetic leads, the conservation of charge current dictates 
that I c ,i = Iq and I Cj % — —Io at the F1|N and F2|N interfaces, which gives 

Mo ~Mo =-(/*o -Mo) = + i;P s v ( 7 ) 

ge 2 
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The free layer is electrically floating, hence there is no net charge current flowing through F3|N interface, I Cj s = 0. 
The spin accumulation in the free layer /if = fi^ — /ij, directed along magnetization m 3 , is governed by the spin 
diffusion equation^ 

which satisfies the following boundary conditions. At the interface the continuity of longitudinal spin current dictates 



V dz 



\ dz J , 



hA 



3 • m 3 



and the vanishing spin current at the end of the ferromagnet implies 



( 

V dz 



dz 



o . 



The solution of eq. © reads 



/*f 0) 



C 3 cosh(^)m 3 .S 



C 3 + a tanh( J- ) cosh( ) 



(9) 



(10) 



(11) 



where £3 = 33(1— P3)/8ir characterizes the contact F3|N and a = hAa^(Ti/(e 2 l s d(cr^+ai)) describes the bulk properties 
of the free layer with arbitrary 1113. The conservation of spin currents in eq. © generates three linear equations that 
determine the spin accumulation S in the normal metal as 



S=n(0,03)(87TlW+W 6 



(12) 



where the vector W& = (0, 2hpIo/e, 0) is the contribution from the bias current and the elements of the symmetric 
matrix 11(3,33) is listed in the Appendix. Equation J3J determines the spin transfer torque acting on the free layer 
magnetization, which can be arranged as 



%53p 

8tt 



(9,93) (8vr 



+ W, 



(13) 



and the components of the matrix T(g,g3) are listed in the Appendix. 



IV. THERMAL STABILITY 



The spin transfer torque rotates the magnetization out of the equilibrium hence increasing the magnetostatic energy 
Ems- The initial magnetization is stable against thermal fluctuations when 

A£ MS > k B T 7 (14) 

where ks is the Boltzmann constant and T the temperature. For an elliptic permalloy film, disregarding any residual 
crystalline anisotropy, the effective field due to the shape anisotropy can be written as 

H e = -fi M s (N x m x , N v my, N z m z ) , (15) 

introducing the saturation magnetization M s and demagnetizing factors N x , N y and N^&. When the magnetization 
is slightly out of plane, the difference between the magnetostatic energy for magnetizations along the hard-axis and 
easy- axis reads AEms = H0VM 2 (N y — N x )/2. For a very flat ellipsoid (the thickness is much smaller than the lateral 
dimensions) and slight ellipticity (large aspect ratio £ w 1), we can expand the demagnetizing factors at £ = 1 such 
that 



(16) 
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where a, b and d are the lengths of easy-axis, hard axis and the thickness of the permalloy film. The aspect ratio is 
defined as £ = b/a. The requirement AEms > kgT gives 

($ 2 + 4£+f)(l-Q 8kB T 

(£ + l) 2 fi M^7T 2 ad 2 ' 1 ' 

The saturation magnetization of permalloy is M s = 8 x 10 5 A m -1 . For thickness d = 5 nm and easy axis a — 200 nm 9 , 
the right hand side of eq. (|17fl is at room temperature approximately e = 8.36 x 10 -3 and therefore stability requires 
that 

£<1-^, (18) 

which suggests that even for almost circular permalloy discs, e.g., £ = 0.9, thermal fluctuations around the equilibrium 
configuration are small. 



V. MAGNETIZATION DYNAMICS 



Here we focus on the free layer magnetization dynamics in the macrospin model. The Landau-Lifshitz-Gilbert 
(LLG) equation modified by the spin transfer torque [eq. 1|13[) ] reads 

1 dm an dm 1 , . 

r- = -m x H e + — m x — + — — L . (19) 

7 dt 7 dt VM S v ' 



Included in the torque term, i.e., Eq. (|13(l . an expression 

<->, = T^jgi 
8%VM S 



r(ff,53) (20) 



appears as an enhancement of the Gilbert damping, 17 : which depends on the direction of the magnetization and shows 
tensor property of the pumping induced damping enhancement. 21 When the conductance at the F3|N contact is much 
larger than the source-drain contacts and the spin flip in the normal metal is negligible, i.e., g 3 3> g and 3> g s f, 
the tensor <r a' converges to2i 



7?iRe 5 J, 



Ti / 1 - m 2 -m x m y -m x m 



4irVM s 

which in the LLG equation reduces to a diagonal matrix 



X llly — lll. X lll, Z 

-m x m y I - ml ~m y m z I , (21) 



nT, x m z —m y m z 1 — mj 



= 1 m £ gl i 

AttVM s 1 ' 

and the coefficient in front of the matrix is exactly the value derived for the single F|N junction 17 . In the same limit, 
the bias-driven term of the torque reads, 

1 -u - ^ ( 7?m* i . ,23, 



VM„ 2VMJe 



-m v m z 



In the following discussions, we denote the enhanced Gilbert damping parameter as a — an + a' , where a' is the 
diagonal entry in ea. (|22|l . The LLG equation then reads 

1 dm a dm 1 , , 

-— = -mxH e + -mx — + — — L b . 24 

■y dt 7 dt V M s 

For ultrathin permalloy films, without external field and crystalline anisotropy, the magnetization is confined in the 
plane by the shape anisotropy field given by eq. (|15|1 . Equation (|24(l is a nonlinear differential equation that can be 
reformulated as 

^=f(m,/ ) (25) 



5 



where f(m, Jo) is a vector function of magnetization m and bias current To- According to the theory of differential 
equations^, we find two "equilibrium points" at which dm/dt vanishes riii = (1, 0, hpIo/[2efioVM^(N z — N x )]) and 
rfi2 = (0,1,0). Expanding eq. 124(1 at point rii2 and keeping only the first-oder derivatives with respect to the 
magnetization, i.e., 



dm _ / Of 
dt \ dm 



where di/dm is a matrix with elements given by dfi/drrij. Equation l|26[l has non-zero solution when 



det 



di 
dm 



mo 



(26) 



(27) 



This determines the critical current that is necessary to obtain the maximum in-plane rotation, i.e., ir/2: 



Ic - n P (28) 

The LLG equation augmented by the spin transfer torque for the present configuration suggests a two-state behavior 
of the magnetization: Below the critical current I c , the magnetization is pushed out of the initial position (easy axis), 
then undergoing damped precessions and finally stops along the easy axis but with a small z-component, i.e., the 
equilibrium given by rhi. At that position, the demagnetizing field is balanced by the spin torque. Above the critical 
current, the magnetization precesses out of the easy axis and rotates to the hard axis without any precession. 

We simulate the magnetization dynamics for a polarization p = 0.4 of the F3|N and a real part of the mixing 
conductance Reg^A^ 1 = 4.1 x 10 15 cm 2 m& The long semi-axis, short semi-axis and the thickness of the permalloy 
island are a = 200 nm, b = 190 nm, and d = 5 nm, respectively. The calculated demagnetizing factors are N y — 0.0224 
and N x = 0.0191iSl The single-spin density of states in the normal metal is vdos — 2.4 x 10 28 eV^ 1 m _3 iS The 
bulk value of the Gilbert damping parameter is ao = 0.006 and the calculated enhancement of Gilbert damping is 
a' = 0.015. 17 According to eq. I|28[) . for the above dimensions, the critical current to achieve <fi — ir/2 is I c — 139 mA, 
which agrees well with the numerical results. Below the critical current, e.g., Iq = 30 mA, the equilibrium z- 
component determined by the expression of rhiis 0.0087, which also agrees with the numerical results shown in Fig. 

The trajectory of the magnetization when suddenly switching on the bias current Iq = 30 mA is depicted in Fig. 
121 The magnetization starts from the easy axis (point I in the figure), undergoes a damped oscillation and finally 
stops at point F, where the spin transfer torque induced by the spin accumulation in the normal metal is balanced 
by the torque generated by the anisotropy field. Fig. [7| shows the trajectory of the magnetization under switching on 
a the bias current Iq = 160 mA, which is above the critical current. Figures and H3 are the time dependence of the y 
and z-components of the free layer magnetization. These figures indicate that the magnetization response to a large 
current is close to a step function. A smaller size of the permalloy film requires a smaller critical current, as indicated 
by eq. I|28|) . In the above simulation, we did not take into account the effect of a finite RC time for switching on the 
bias current. A longer rising time of the bias current implies that it takes longer before the magnetization reaches the 
steady state position. But the magnitude of the critical current does not depend on how the bias current transient. 

We finally note that with the dimensions chosen here, the bias currents generate a significant in-plane 0rsted field 
that may interfere with the spin-torque effect. It can be avoided, e.g., by spatially separating the free layer from the 
current path (but within the spin-flip diffusion)^ or by generating a neutralizing 0rsted field by a neighboring circuit 
(suggested by Siegmann). 

The advantage of the proposed device mainly comes from the two-state behavior separated by the critical current, 
which can be utilized as the and 1 states in current controlled memory elements. We notice that after the magneti- 
zation being switched to the hard axis, only small current is needed to maintain the position stable against thermal 
fluctuations. Another possible application could be the implementation of such device into spin-torque transistors^ 
to achieve the gain of current since the angle of the magnetization in the above device is tunable by the bias. The 
magnetization can be also used as a spin battery that is "charged" in the high energy state (hard axis) and relaxes 
a spin current into the normal metal when relaxing to the ground state (easy axis). The induced spin accumulation 
then creates voltage difference over the source and drain contacts. 



VI. CONCLUSIONS 



In this article, the magnetization dynamics of a spin-flip transistor has been studies in the macrospin LLG equation 
combined with MECT. We found a two-state behavior of the free layer magnetization controlled by the current 
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FIG. 2: The trajectory of the magnetization with the bias current lo = 30 mA, which is below the critical current. The 
magnetization initially aligned along easy axis (rr-axis) and after the damped oscillation it stops along the easy axis with small 
out-of-plane component. 
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FIG. 3: The x-component of the magnetization vs time (in ns). The bias current is 30 mA. 
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FIG. 4: The 2-component of the magnetization vs time (in ns). The bias current is 30 mA. 
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FIG. 5: The y-component of the magnetization vs time. The bias current is 160 mA, which is above the critical current. 
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FIG. 6: The z-component of the magnetization vs time, with bias current 160 mA that is above the critical current. 




0.025 




FIG. 7: The trajectory of the magnetization under the bias above the critical current, 7o = 160 mA. The magnetization 
initially aligned along easy axis (point I). 
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induced spin transfer torque and spin pumping. The two regimes are separated by a critical current, below which 
the magnetization undergoes a damped oscillation and stops along the easy axis with small z-component. Above the 
critical current, the magnetization rotates to the hard (y-) axis without precession. The critical current is found to 
depend on the size of the free layer, the aspect ratio of the ellipsoid, and the source-drain contact polarizations. The 
thermal instability analysis indicates that at room temperature the predicted effects are visible even for very large 
aspect ratios. 
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APPENDIX A: SPIN ACCUMULATION AND SPIN TRANSFER TORQUE 



The elements of the symmetric matrix tl(g, g 3 ) in eq. (|12f) are listed in the following 

Hii =[V393 + 2g si + 2(1 - p 2 ) 5 ][2 w + 77333 + 2gsf - <? 3 (% - (1 - P§)(1 - A 3 ))m=]/g (Al) 

- (2775 + r/3.93 + 2. 9sf ) 53 b73 - (1 -p§)(l - &3)}m 2 y /g 

U 12 =n 21 = (2 V g + 773.93 + 2g Bi )g 3 [r] 3 - (1 - - A 3 )]m x m y /g (A2) 

n 13 =n 3 i = [77333 + 2 3sf + 2(1 - p*)g]g 3 [ m - (l - pi) (1- A 3 )]m x m z /g (A3) 

n 22 =(277.g + 773.93 + 2g sf )[277.9 + 77353 + 2g si - .93(773 - (1 - p§)(l - A 3 ))(m£ + m*)]/0 (A4) 

n 23 =n 32 = (277.9 + 773.93 + 2g si )g 3 [r ]3 - (1 - pg)(l - A 3 )}m y m z /g (A5) 

n 33 =[r? 3 <?3 + 2 5sf + 2(1 - p^g^g + m g 3 + 2g sf - g 3 (r] 3 - (1 - p§)(l - A 3 ))m 2 x }/g (A6) 

- (2t7 5 + V393 + 2g si )g 3 [m - (1 - - A 3 )]m 2 y /g 



where we have introduced the following notation 

g = (2to + 773.93 + 2 9s f)[(7 M3 + 2g- 2p 2 g + 2g st )(2r)g + 2g sl + (1 - - A 3 ) 53 ) (A7) 
+ 2(p 2 - 1 + 77)3(773 - (1 - _ A 3 )) 53 m2] , 

and 

A3 = C 3 + ^ta C nh(d// sd ) ■ (A8) 
The matrix contained in the expression of spin transfer torque eq. (|13Jl has the following components 

Tn = [53.93 + 2g ai + 2(1 - p 2 )g][2 V g + 77393 + 2.g sf - 53(773 - (1 -p\){\ - A 3 ))m 2 z }/g (A9) 

- (2775 + 77353 + 2ftrf)[(»fe ~ (1 - A 3 ))g 3 m 2 y + ( m g 3 + 2g s{ + 2(1 - p 2 )g)m 2 x ]/g 

r u = -(255 + 53.93 + 2g ai )[2r)g + 2g sf + (1 - p§)(l - A 3 )fl3]m x m,,/e (A10) 

Ti 3 - -[53.93 + 2.9sf + 2(1 - p 2 )g][2 V g + 2g sf + (1 - p§)(l - A 3 )g 3 ]m x m z /g (All) 

T 2 i = -(255 + 5353 + 2g ai )[2g B{ + 2(1 - p 2 )g + (1 - p 2 )(l - A 3 )g 3 }m x m y /g (A12) 

T 22 = (255 + 5353 + 25sf)[25.9 + 2 3 sf + (1 - P§)(1 - A 3 )g 3 ]{l - rr^)/g (A13) 

T 23 = -(255 + 5353 + 2 5sf )[2 ffsf + 2(1 - p 2 )g + (1 - p§)(l - A 3 )g 3 ]m y m z /g (A14) 

T 3 i = -[53.93 + 2g st + 2(1 - p 2 )g][2i lg + 2g si + (1 - p 2 )(l - A 3 )g 3 ]m x m z /g (A15) 

T 32 = -(255 + 5353 + 2 5sf )[255 + 2 3sf + (1 - p 2 3 ){l - A 3 )g 3 }m y m z /g (A16) 

T 33 = [5353 + 2 5sf + 2(1 - p 2 )g] [255 + 5353 + 2 3s f - 53(53 - (1 - p|)(l - A 3 ))m 2 x ]/g (A17) 

- (255 + 53.93 + 2.9 sf )[(53 - (1 - p§)(l - A 3 ))g 3 m 2 y + (53.93 + 2 5sf + 2(1 - p 2 )g)m 2 z ]/g . 



